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THE  STATISTICAL  ANALYSIS  OF  THE  DISTRIBUTION  OF  SNOW  COVER  DEPTH  IN  SMALL 

MOUNTAINOUS  AREAS 

By:  A.  A.  Chirkova 

Questions  of  Improving  the  efficiency  of  the  snow-measuring  system,  despite 
significant  experience  recently  accumulated  (I,  4,  5,  7,  14,  15),  still  remain 
pressing  ones:  first,  because  of  the  Increase  In  Information  about  the  statis- 
tical structure  of  the  snow  reserve  field,  and  second,  because  of  the  necessity 
of  constructing  a system  using  new  means  of  data  transmission,  and  finally.  In 
connection  with  the  practical  requirements  of  differentiated  estimates  of  snow 
reserves . 

This  work  makes  an  attempt  to  estimate  the  amount  of  snow  cover  In  areas 
comparable  to  the  dimensions  of  a snow  (observation)  post  using  statistical 
methods.  The  original  material  wets  the  results  of  special  snow-meatsurlng  surveys 
conducted  in  1967,  1971  and  1972  in  certain  basins  of  western  Tyan'-Shan'. 

The  principles  of  choosing  the  area  consisted  In  an  attempt  to  find  the 
most  level  parts  of  a slope  and  to  hold  the  effect  of  characteristics  of  snow 
accumulation  In  any  particular  basin  to  a minimum.  Measurements  of  the  depth 
of  snow  cover  were  made  along  lines  50  m long  rigidly  held  in  place:  in  6 
C£ises  edong  four  mutually  Intersecting  lines,  and  In  4 cases  along  the  perimeter 
of  a square.  Distance  was  fixed  with  the  aid  of  a measuring  tape  every  20  or 
50  cm.  A total  of  28  series  of  snow  surveys  was  carried  out. 

Statistical  analysis  of  the  obtained  data  started  with  a check  of  the 
assumption  of  homogeneity  of  the  studied  field. 

A field  of  mean  values  of  the  height  of  the  snow  cover,  followed  by  the 
dispersion  field,  was  primarily  Investigated.  The  Student  criterion  was  used 
to  test  the  equality  (homogeneity)  of  two  means 


/I,  + /I, 


where  X2^  and  X2  - the  random  means;  and  02  - random  dispersions;  n^^  and  n2  - 
number  of  memebers  of  the  first  and  second  samples;  (n^  + n2  -2)  - number  of 
degrees  of  freedom. 

For  normally  distributed  samples  at  the  q il-level  of  significance,  deviation 
between  the  centers  of  distribution  is  considered  significant  if  the  following 
condition  Is  satisfied 


> t 


The  boundaries  of  the  critical  region  have  been  tabulated  and  cited  In  works 
(8,  13)  for  a Student  distribution. 

The  results  of  the  test  given  In  Table  1 show  that  In  most  cases  the  random 


student  statistic  exceeds  the  critical  value  at  the  5->percentile  level  of  sig- 
nificance (tg^  = 1.96). 

We  note  that  the  use  of  the  Student  criterion  is  valid  with  the  condition 
of  equality  of  dispersions  and  can  lead  to  erroneous  conclusions  when  this 
requirement  is  not  observed.  Therefore,  in  addition  to  testing  the  equality 
of  dispersions,  the  affiliation  of  the  means  to  a single  set  was  randomly 
tested  again  according  to  the  criteria  of  signs  not  limited  by  the  equality 
of  dispersions.  As  in  the  preceding  case,  the  results  of  the  test  validated 
the  hypothesis  of  homogeneity  of  mean  values. 


The  F-statistic  of  Fischer  was  used  to  test  the  equality  of  dispersions. 
This  statistic  makes  it  possible  to  determine  the  significance  of  deviations 
among  the  cuialyzed  dispersions.  The  difference  is  recognized  as  significant 
if  the  following  condition  is  satisfied 


F > F . 
q 


The  determination  of  the  criterion  which  depends  solely  on  the  number  of 
degrees  of  freedom  pertains  to  calculating  the  value 

®1 

F = 


2 

where  s|  is  always  the  greater  of  the  examined  values.  In  this  case  the 
estimate  of  dispersion  is  determined  for  each  measured  line  and  for  the  entire 
realization  as  a whole.  The  5-percentile  level  is  taken  as  the  calculation 
level . 


The  calculated  values  of  the  Fischer  criterion  for  most  cases  were  over 
the  critical  ones,  which  indicates  the  existence  of  differences  of  dispersions 
(Table  1).  It  is  physically  difficult  to  explain  the  obtained  results,  inas- 
much as  the  conditions  of  snowfall  and  snow  deposit  within  the  confines  of  the 
experimental  areas  are  practically  unaltered.  Analysis  of  the  Fischer  statsitics 
obtained  for  each  measuring  line  detects  an  effect  of  a certain  factor.  To  our 
view,  this  could  be  the  character  of  the  underlying  surface.  Because  of  this, 
the  significance  of  its  effect  on  the  random  nature  of  the  sample  was  investi- 
gated. The  estimate  was  made  with  the  aid  of  the  successive  differentials 
methods  (the  Neiman  relationship),  in  which  the  random  dispersion  is  estimated 
in  two  ways: 

s*=-^E  (jc,— T)*, 

" (-1 


! 

“ 2(n-l) 


S (-*<+1  •*!)*• 


2 2 

Estimates  c and  s react  differently  to  a change  in  the  center  of  the  distri- 
bution. The  former  is  less  sensitive  to  the  trend  of  the  mean  value,  while  the 
latter  sharply  increases  in  the  presence  of  even  a weak  trend. 
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Key  for  Table  1: 


1 - Table  1.  Homogeneity  test  of  average  values  of  snow  depth  with  the  aid  of 

the  Student  criterion  (1)  and  of  dispersions  with  the  aid  of  the  Fischer 
criterion  (2)  at  the  5-peroentile  significance  level. 

2 - date  of  measurements 

3 - measuring  line 

4 - Kyzylcha 

5 - Naugarzan 

6 - Oygaing 

7 - Lednik  Severtsova 

(Note:  commas  should  be  read  as  decimals.) 
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When  the  sample  has  a large  volume  i.he  relationship 


2 

s 


approximates  the  normal  and  therefore,  one  can  use  the  quantiles  of  normal  dis- 
tribution, whose  tables  are  given  in  a number  of  works,  to  determine  the  boundaries 
of  the  critical  region 


The  trend  is  significant  if  t < tg  at  the  assigned  level  of  significance. 

In  the  examined  case,  when  q = and  n = 100,  the  critical  value  t = 0.835, 
and,  as  is  evident  from  Table  2,  it  significantly  exceeds  the  calculated  values 
which  are  near  zero  everywhere,  which  is  a characteristics  of  the  correlated 
initial  data  (12).  Consequently,  the  hypothesis  of  the  random  nature  of  the 
original  field  contradicts  the  experimental  data  and  should  be  discarded.  Un- 
doubtedly, this  is  the  reason  why  the  homogeneity  hypothesis  was  overturned, 
inasmuch  as  it  is  known  that  the  hypothesis  does  not  take  into  account  the 
intraset  relations  and  leads  to  taking  patently  homogeneous  data  as  hetero- 
geneous ones  (12).  Therefore,  subsequently,  one  should  designate  the  interval 
between  measurements  in  order  to  ensure  independence  of  the  original  data. 

Some  of  the  estimates  used  above  (the  parameteric  criteria)  depends  sig- 
nificantly on  the  normality  of  the  distribution  and  for  this  reason  require 
a basis  for  the  validity  of  their  use.  At  the  same  time,  opinions  which  exist 
in  the  literature  concerning  the  form  of  snow  cover  distribution  are  extremely 
contradictory.  Some  authors  (2,  3,  11)  assume  that  the  snow  has  a normal,  or 
at  least  normal  distribution,  while  others  (10,  16)  include  the  snow  cover 
among  phenomena  with  a clearly  pronounced  asymmetry. 

The  task  of  the  investigation  the  authors  conducted,  first,  was  to 
establish  the  law  of  distribution  in  small  areas,  and  second,  to  analyze  the 
existence  of  a deviation  of  the  original  theories  from  the  normal  law.  In 
order  to  describe  the  type  of  distribution  of  the  snow  cover,  the  widely 
popular  Pearson  curves  were  used.  For  such  curves,  methods  and  systems  of 
calculation  are  well  developed.  In  our  case,  smoothing  of  the  empirical 
frequencies  was  carried  out  according  to  the  system  presented  in  (8). 

The  calculations  were  carried  out  on  a computer  according  to  a program 
compiled  by  G.  Ye.  Glazyrin  (SARNIGMI),  which  provides  not  only  for  the  choice 
of  the  type  of  Pearson  curve,  but  also  for  a comparison  of  the  empirical  fre- 
quencies with  those  chosen  with  the  aid  of  the  Pearson  agreement  criterion. 

The  program  also  provides  for  analog  operations  with  a normal  curve  additionally 
entered  into  the  auialysis.  In  parallel  with  the  Pearson  criterion,  the  hypo- 
thesis of  the  correspondence  of  empirical  data  to  the  suggested  theoretical 
law  was  tested  according  to  the  Kolmogorov  criterion,  since  the  Pearson 
criterion  significantly  differs  from  the  quite  arbitrarily  accomplished  grouping, 
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and  furthermore,  the  use  of  the  maximum  law  of  distribution  with  small  volumes 
of  the  sample  will  lead  to  the  recognition  of  a zero  hypothesis  even  for  sig- 
nificantly different  distributions,  while  for  large  samples,  on  the  other  hand, 
it  frequently  refutes  the  zero  hypothesis  (6). 

In  both  cases,  the  assumption  was  made  for  the  zero  hypothesis  that  the 
examined  sample  is  distributed  according  to  the  hypothetically  theoretical 
law. 

The  following  expression  is  used  as  a measure  of  deviations  in  the  Pearson 
criterion 

i 

where  nj^  - empirical  frequencies;  rii  - smoothed  frequencies;  1 - number  of 
gradations. 

At  the  5/l-level  of  significance,  the  critical  region  for  testing  the 
hypothesis  has  the  following  form: 

The  results  of  the  check  are  given  in  Table  3.  In  analyzing  the  tabular 
data,  it  is  impossible  to  give  preference  to  either  asymmetrical  or  symmetrical 
distribution:  the  appearance  of  both  the  first  amd  second  curves  is  practicably 
equiprobable.  Thus,  of  50  cases,  Pearson  curves  of  the  first,  third,  fourth 
and  sixth  types  (asymmetrical)  are  encountered  22  times,  amd  in  28  cases  the 
criterion  satisfies  the  requirements  of  curves  of  the  symmetrical  type  (2 
and  7).  Both  kinds  of  distribution  are  observed  at  all  experimental  sites, 
independent  of  the  uniqueness  of  these  sites.  The  asymmetrical  distribution 
predominates  only  in  regions  with  pronounced  wind  activity.  One  also  notes 
a tendency  of  more  frequent  mainifestation  of  the  asymmetrical  curves  in  the 
snow  thaw  period. 

The  Kolmogorov  criterion  is  based  on  a comparison  of  the  integral  distri- 
bution functions.  The  measure  of  deviation  here  is  the  greatest  difference 
between  the  empirical  and  theoretical  functions  relative  to  absolute  magnitude. 


D = max  P!  - P.  . 

• 1 1 ' 

The  statistic  \ = D-^nc"  with  a critical  probability  of  K(x)  = 0.5  means  that 
deviations  between  the  examined  distributions  are  due  to  random  factors  and 
are  not  significant. 

The  results  of  the  test  are  given  in  Table  3.  According  to  this  criterion, 
the  symmetricad  Pearson  curves  are  encountered  41  times  and  the  asymmetrical  ones 
are  encountered  36  times.  It  is  more  important  to  note  that  if  a normal  distri- 
bution was  found  in  44  samples  of  the  103  investigated  ones  according  to  the 
Pearson  criterion,  then  during  the  use  of  the  Kolmogorov  criterion  the  model  of 
normal  distribution  can  be  accepted  in  78  cases. 
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Key  for  Table  2: 

1 - Table  2.  A check  for  the  random  nature  of  snow  cover  depth  measurements 

by  the  successive  differential  method. 

2 - date  of  measurements 

3 - measuring  line 

4 - Kyzylcha 

5 - Naugarzan 

6 - Oygaing 

7 - Lednik  Severtsova 

8 - Dukant-1 

9 - Dukant-2 
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This  conclusion  is  a very  important  one  in  a practical  regard,  since  it 
imparts  a justification  to  the  calculations  done  above  and  makes  it  possible 
to  determine  what  volume  of  observations  one  should  have  to  find  the  mean  values 
with  an  assigned  level  of  error  and  probability. 

2 

In  comparing  the  results  obtained  by  means  of  calculating  the  x criterion 
and  the  Kolmogorov  criterion  for  the  hypothesis  of  correspondence  of  the  sample 
data  to  the  theoretical  law,  it  is  easy  to  see  that  on  the  whole  both  criteria 
agree  well.  It  is  only  for  samples  with  the  lowest  probability  = 0.30) 

of  statistic  X that  the  criterion  produces  a negative  result. 


We  now  proceed  to  determine  the  minimum  number  of  measurements  according 
to  a formula  known  in  statistics: 


It  — 


where  tp  is  a value  determined  at  an  assigned  probability  P in  tables  of  values 
of  the  probability  integral;  Cy  is  a measure  of  changeability  (the  coefficient 
of  changeability);  e is  the  permissible  relative  error.  Usually,  in  the  practice 
of  scientific  investigations,  the  probability  value  is  assumed  to  be  P = 0.95 
or  P = 0.99,  while  the  permissible  measurement  error  is  established  depending 
upon  the  nature  of  the  phenomenon.  As  the  investigations  of  I.  A.  Il’in  (SARNIQII) 
showed,  in  90$  of  the  cases  the  possible  error  in  determining  the  average  depth 
of  the  snow  does  not  exceed  10$  in  the  period  of  build-up  of  the  snow  cover  with 
a certainty  value  of  0.95.  In  March  the  error  increases  and  in  the  same  90$  of 
the  cases  already  comprises  about  20$.  These  investigations  were  carried  out 
according  to  the  data  of  snow -measuring  posts  in  mountain  basins  of  Central  Asia. 
The  results  of  the  calculation  at  a certainty  level  of  0.95  in  an  error  of  0.5 
in  one  case  and  0.10  in  another  are  given  in  Table  4.  It  is  entirely  obvious 
that  one  cannot  have  a uniform  solution  when  one  selects  the  error  of  determining 
the  mean.  With  precisely  the  same  number  of  measurements,  the  error  of  averaging 
will  first  be  different  at  different  times  of  snow  accumulation,  and  second,  will 
be  different  for  areas  with  a different  degree  of  terrain  brokenness : the  greater 
the  degree  of  brokenness  of  the  terrain,  the  more  measurements  should  be  made  in 
the  area  in  order  to  obtain  an  average  with  an  assigned  level  of  accuracy. 

Figure  1 visibly  confirms  the  above.  Here,  on  the  axis  of  the  ordinate,  one 
has  the  necessary  number  of  measurements  calculated  when  e = 0.10  and  P = 0.95 
for  the  period  of  build-up  of  the  snow  cover;  on  the  axis  of  the  abscissa  one 
has  the  coefficient  of  changeability  of  terrain  elevations  over  an  arbitrarily 
selected  area.  The  latter  are  values  accepted  as  the  index  of  terrain  broken- 
ness measured  in  the  field.  And  there  is  still  another  important  detail.  In 
errors  with  wind  redistribution  of  the  snow  cover  (for  example,  in  an  area 
located  in  the  Naugarzan  River  basin),  the  effect  of  the  underlying  surface 
is  smoothed  by  changeability  created  by  wind  snow  transport.  In  this  case  the 
number  of  necessary  measurements  in  such  areas  increases  many  times  in  comparison 
with  areas  which  are  not  subjected  to  the  wind  effect. 
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Key  for  Table  3: 


1 - Table  3.  Testing  the  hypothesis  of  distribution  of  measurements  of  the  snow 

cover  depth  with  the  aid  of  the  (D  and  Kolmogorov  (2)  criteria. 

2 - date 

3 - number  of  measuring  line 

4 - Pearson  curve 

5 - number  of  measurements 

6 - type  of  curve 

7 X^sam 

8 - estimate  Hq 

9 - normal  curve 


10 

- 

X^sam 

11 

- 

estimate 

Hq 

12 

- 

Pearson 

curves 

13 

- 

normal  curve 

14 

- 

estimate 

Ho 

15 

- 

estimate 

Ho 

16 

- 

Dukant-1 

17 

- 

Dukant-2 

18 

- 

yes 

19 

- 

no 

In  suimarlzlng  the  results  of  the  Investigation,  one  can  assert  the 
following:  when  one  observes  the  interval  of  independence  between  meaisure- 
ments,  in  which  the  intraset  correlation  is  held  to  a minimum,  the  distri- 
bution of  snow  cover  at  the  snow-measuring  post  should  be  assumed  to  be  homo- 
geneous with  a nearly  normal  distribution.  In  this  case,  in  the  most  general 
form  (the  area  is  either  horizontal  or  has  a slight  slope  and  is  even),  when 
one  has  a relative  error  of  the  average  of  0.10  it  is  sufficient  to  make  5-10 
measurements  at  the  snow-measuring  post  in  order  to  obtain  a mean  value  with  a 
probability  of  0.95,  but  2dready  when  one  has  a relative  error  of  0.05,  the 
number  of  observations  should  be  increased  to  10  - 15. 


TafiAiiua  4 
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Key  for  Table  4: 


1 - Table  4.  The  minimum  number  of  measurements  necessary  to  obtain  averaige 

values  of  snow  depth  with  a probability  of  0.95. 

2 - date  of  measurements 

3 - measuring  line 

4 - tot^d 

5 - Dukant-1 

6 - Dukant-2 

7 - Kyzylcha 

8 > Naugarazan 

9 - Oygaing 

10  - Lednik  Severtsova 


Figure  1.  The  relationship  of  the  number  of 
measurements  with  the  brokenness  terrain. 
Areas  in  basins:  1 - Kyzylcha  River;  2 - 
Oygaing  River;  3 - Kashkadar'ya  River;  4 - 
Dukant  River. 
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i 

I 

I 

I In  selecting  the  area  of  the  snow-measuring  post,  it  is  necessary  to  avoid 

I places  with  wind  snow  redistribution  (drifts,  denuded  spots),  and  if  necessary 

to  select  an  even  area  taking  into  account  the  nature  of  terrain  brokenness.  ^ 
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